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ABSTRACT 


This  report  presents  survey  material  on  the  effects  of  the  reentry  environment  on  the  voltage 
breakdown  characteristics  of  antennas.  Most  of  the  important  effects  are  fairly  well  under¬ 
stood,  in  that  theoretical  models  exist  and  experimental  data  are  available  in  sufficient 
quantity  to  allow  predictions  to  be  made  with  confidence.  However,  the  effect  of  high  gas 
temperature  is  not  sufficiently  understood,  although  recent  experimental  data  have  indicated 
a  mechanism  which  may  prove  to  be  important  in  controlling  breakdown.  This  mechanism 
is  strongly  dependent  on  vibrational  excitation  and  dissociation  of  the  gas  molecules. 
Predicting  the  vibrational  excitation  state  in  a  hypersonic  flow  field  is  not  presently  possible 
in  a  detailed  microscopic  analysis,  but  an  approximate  solution  is  possible. 
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SECTION  1 
INTRODUCTION 

The  study  of  which  this  report  is  a  part  is  concerned  with  the  effects  of  the  reentry  environ¬ 
ment  of  a  slender  hypersonic  vehicle  on  the  transmission  of  ECM  signals.  The  most 
important  aspect  of  this  study  is  the  effects  of  reentry  on  antenna  breakdown.  In  response 
to  technical  direction,  the  effort  during  this  reporting  period  was  concentrated  on  a  survey 
of  reentry  effects  on  antenna  breakdown  and  included  visits  to  the  facilities  of  some  of  the 
workers  in  the  field. 

Part  of  the  scope  of  effort  involves  fluid  mechanics.  It  turns  out  that  one  of  the  most 
important,  but  least  easily  analyzed,  effects  of  reentry  on  antenna  breakdown  is  the  effect 
of  high  air  temperature  in  terms  of  vibrationally  excited  states  of  molecular  constituents. 

As  a  result,  and  in  response  to  direction,  the  fluid  mechanics  efforts  were  devoted  to  the 
initial  phases  of  a  review  of  vibrational  excitation  effects  in  a  nonequilibrium  flow  field 
situation. 


SECTION  2 

FLOW  FIELD  ANALYSIS 

In  accordance  with  the  recommendations  which  were  made  subsequent  to  the  first  technical 
direction  meeting,  the  flow  field  literature  survey  has  been  rescheduled  to  be  completed  by 
the  end  of  the  third  contract  quarter.  Following  the  suggestions  which  were  made  at  the 
meeting,  it  is  intended  that  during  this  survey,  special  emphasis  will  be  directed  toward 
vibrational  nonequilibrium  effects  on  the  flow  field  properties.  Part  of  this  literature  survey 
has  been  completed  and  is  reported  below.  Also  reported  below  is  the  development  of  a 
relatively  simple  method  which  will  be  used  to  determine  the  extent  of  the  effects  of 

vibrational  nonequilibrium  on  the  flow  field  calculations  which  were  presented  in  the  first 
quarterly  report. 

2«1  FLOW  FIELD  LITERATURE  SURVEY 

During  the  course  of  the  flow  field  survey  a  report  entitled  'Interaction  of  Chemical  Reaction 
With  Vibration  or  Translation  in  Nonequilibrium  Processes"  by  E.  Schwartz  of  GE/RSD36* 

was  reviewed  and  found  to  pertain  directly  to  this  survey.  Most  of  the  following  material 
has  been  taken  from  this  report# 

When  air  passes  through  a  shock  the  Maxwellian  velocity  distribution,  disrupted  by  the 
highly  nonequilibrium  conditions  in  the  shock,  is  restored  within  the  time  of  a  very  few 
collisions  ,  and  the  translational  degrees  of  freedom  of  the  gas  are  equilibrated  almost 
instantaneously  at  the  new,  higher  temperature. 

The  relaxation  time  for  rotation  is  somewhat  longer,  about  a  dozen  collisions2  being  required 
to  bring  the  molecular  rotation  into  equilibrium,  i.  e. ,  to  establish  the  Maxwell-Boltzmann 
distribution  of  rotational  energy  levels  characteristic  of  the  new  translational  temperation. 

The  time  required  for  this  is  still  small  enough  to  justify  its  neglect. 


♦References  appear  in  Section  6. 
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The  relaxation  of  vibration  is  much  slower;  for  instance,  the  approximate  number  of  > 

collisions  required  to  establish  vibrational  equilibrium  for  Ng  and  02  at  5000°C  is  3000  and 
200,  respectively3.  Vibrational  disequilibrium  must  be  taken  into  account  in  reducing  shock 
tube  data  and  in  making  accurate  calculations  cf  hypersonic  flows. 

Lastly,  the  slowest  processes  in  reactive  flows  are  the  chemical  reactions  themselves.  In 

earlier  work  these  rates  were  assumed  to  be  infinite;  i.  e. ,  the  concentrations  of  various 

species  were  obtained  from  equilibrium  constants.  But  for  high-speed  flows  this  assumption 

is  entirely  inadequate.  For  instance,  species  concentrations  calculated  on  the  assumption  of 

4 

chemical  equilibrium  may  be  orders  of  magnitude  in  error  .  Thus,  the  concentrations  of 
the  various  species  must  be  found  from  a  set  of  differential  equations,  one  for  each  reaction 
considered.  Because  the  Lav'  of  Equipartition  of  Energy  is  not  applicable  in  nonequilibrium 
processes,  the  various  energy  modes  are  strongly  coupled  by  the  rates  at  wnich  energy  is 
stored  or  transferred.  The  rate  "constants"  which  are  used  in  the  above-mentioned 
differential  equations  are  particularly  sensitive  to  these  coupling  effects. 

In  all  the  studies  reviewed  below,  effects  of  viscosity,  heat  conduction,  and  diffusion  have 
been  ignored.  Except  in  the  boundary  layer  and  in  the  interior  of  the  shock  wave,  the 
various  gradients  are  small  enough  to  justify  this  neglect  of  transport  phenomena. 

2. 1. 1  PERTURBATION  OF  MAXWELLIAN  VELOCITY  DISTRIBUTION  BY  CHEMICAL 
REACTION 

It  has  long  been  recognized  that  any  chemical  reaction  taking  place  at  a  finite  rate  must 
necessarily  disturb  the  Maxwellian  distribution  of  translational  velocity.  Since  the  reaction 
is  brought  about  by  the  more  energetic  collisions,  the  reaction  may  deplete  the  high-energy 
tail  of  the  velocity  distribution.  On  the  other  hand,  for  reactions  with  high  activation  energy 
and  low  steric  factors  (low  probability  for  reaction  even  for  sufficiently  energetic  collisions!, 
most  collisions  will  not  be  reactive  but  wfll  help  maintain  translational  equilibrium,  and  the 
departure  from  the  Maxwellian  distribution  should  be  negligible.  A  number  of  investigators,  *> 

by  making  perturbation  solutions  of  a  suitably  generalized  Boltzmann  equation,  have  made 
quantitative  estimates  of  this  interaction  of  chemical  reaction  and  translational  disequilibrium.  f 
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Prigogene  and  Xhrouet  consider  the  initial  rate  of  the  reaction  A  +  A-^B  +  C  for  molecules 

of  spherical  symmetry  and  no  internal  degrees  of  freedom.  The  influence  of  B  and  C  is 

ignored  and  the  heat  of  reaction  disregarded.  It  was  also  postulated  that  the  cross  section 

6 

for  reaction  was  small.  The  method  of  Chapman  and  Enskog  was  used  for  solving  a 
modified  Boltzmann  equation  with  the  velocity  distribution  taken  as  f  =  f  (1  +  01  where  fQ  is 
the  equilibrium  (Maxwellian)  distribution  and  0,  a  perturbation.  It  was  found  that  both  0  and 
the  ratio  of  the  reaction  rate  to  its  value  at  equilibrium  depend  only  on  e  */kT,  where  e*  - 
activation  energy.  The  quantitative  results  were  very  sensitive  to  the  choice  of  the  reaction 
cross  section  or,  a  function  of  the  relative  molecular  velocity;  but  the  general  conclusion  was 
reached  that  when  e  */kT  >  5  the  perturbation  of  the  velocity  distribution  and  the  deviation 
from  equilibrium  reaction  rate  are  quite  small  and  that  the  use  of  equilibrium  theory  in 
calculation  of  reaction  rates  was  justified  for  reactions  of  the  type  considered. 

9 

According  to  References  7  and  8,  the  work  of  Curtiss  is  more  general  and  his  approach 
somewhat  different,  but  his  results  are  in  essential  agreement  with  those  of  Prigogene  and 
Xhrouet. 


Prigogene  and  Mahieu^  extend  the  work  of  Reference  5  by  including  the  influence  of  heat  of 

reaction  r  and  find  that  it  can  perturb  the  Maxwellian  distribution  to  an  appreciable  extent. 

They  calculate  the  initial  rate  of  the  reaction  A  +  Br-*-A,  +  B,  where  A  and  A  differ  only 

o  1  o  I 

in  electronic  energy,  and  obtain  the  equation  (for  a  particular  choice  of  reaction  cross 
section): 


v/v  =  1  +  1.2  x.  /€*  (2. 1.1-1) 

o  A  B 

o  o 

where 

P 

v  =  reaction  rate 

f  v  =  equilibrium  v  (i.  e. ,  under  assumption  of  Maxwellian  distribution) 

o 

X  =  mole  fraction 

% 
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This  shows  that  the  reaction  rate  will  be  increased  for  exothermic  reactions  and  that  the 

order  of  the  reaction  is  increased  because  of  the  presence  of  the  term  Xg 

o 
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Present  gives  a  simplified  derivation  of  the  results  of  Prigogene  and  Xhrouet  using  a  more 
realistic  reaction  cross  section.  He  finds  that  even  for  e  */kT  as  small  as  5,  the  reaction 
rate  is  only  8  percent  less  than  that  predicted  from  equilibrium  theory. 

Mahan**  points  out  that,  since  perturbation  of  the  Maxwellian  velocity  distribution  becomes 
important  only  for  reactions  of  low  activation  energy  at  high  temperatures,  the  most  serious 
deviations  from  equilibrium  reaction  rate  may  be  expected  to  occur  in  free  radical  combina¬ 
tion  reactions,  where  little  or  no  activation  energy  is  involved.  He  derives  a  collision  cross 

section  for  radical  combination  reactions  and  then  calculates  0  and  v/v  in  the  manner  of 

o 

Prigogene.  Mahan  concludes  that  the  equilibrium  assumption  is  justified  in  photochemical 
systems  where  free  radical  concentrations  are  low  but  may  not  be  valid  in  flames  and 
detonations  where  radical  concentrations  may  exceed  0. 1  mole  fraction. 

12 

Ross  and  Mazur  present  a  formal  perturbation  solution  of  a  Boltzmann-type  equation  for  a 

bimolecular  reaction  (A  +  B  *  C  +  D),  in  which  the  backward  as  well  as  the  forward  reaction 

is  considered  and  the  molecules  may  have  internal  degrees  of  freedom.  They  conclude  that 

the  usually  assumed  phenomenological  rate  expressions  for  chemical  reactions  (for  example, 

-dn  /dt  =  k„  n  n  -knn,  for  the  reaction  A  +  B  ~  C  +  D)  are  strictly  valid  only  in  the  lowest 
a  fab  red 

order  of  the  perturbation.  Higher  order  perturbations  show  the  rate  "constants"  to  be  a 
function  of  affinity  and  hence,  implicitly,  of  the  time. 

2. 1.2  PERTURBATION  OF  MAXWELL-BOLTZMANN  DISTRIBUTION  OF  VIBRATIONAL 
LEVELS  BY  DISSOCIATION 

Since  a  molecule  in  a  state  of  high  vibrptional  excitation  is  more  likely  to  undergo  dissociation 
on  collision  than  one  of  lesser  excitation  (with  the  same  translational  energy),  a  dissociation 
reaction  (A  +  M~A  +  A  +  M)  will  disturb  the  Maxwell-Boltzmann  (equilibrium)  population 

dt 

of  the  vibrational  levels;  and  the  reaction  rate  will  be  smaller  than  that  calculated  on  the 
basis  of  vibrational  equilibrium. 
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In  order  to  study  vibrational  level  populations  and  dissociation  rates  under  nonequilibrium 
conditionst  the  investigators  cited  below  consider  a  dilute  solution  of  oscillators  (diatomic 
molecules  A0)  in  an  inert  gas  M  which  actf*  as  a  constant  temperature  heat  bath  (i.  e. » 

Li 

effects  discussed  in  Section  2. 1. 1  need  not  be  considered).  The  concentration  of  oscillators 
Is  so  small  that  only  collisions  with  the  heat  bath  molecule  need  be  taken  into  account.  This 
leads  to  a  set  of  N  simultaneous  linear  rate  equations  in  the  Xj  (t)  (N  =  number  of  vibrational 
levels;  X|  =  fraction  of  oscillators  in  level  i),  whose  coefficients  are  the  probabilities  per 
unit  time  of  transition  from  one  level  to  another  or  of  dissociation  from  that  level.  This  set 
of  differential  equations  is  known  variously  as  the  "master  equation"  or  the  "transport 
equations.  "  Methods  are  available  for  their  solution  which  do  not  require  that  the  system 
be  close  to  equilibrium,  as  do  the  perturbation  methods  of  Section  2. 1. 1.  It  should  p 
pointed  out,  however,  that  generally  these  methods  involve  finding  the  eigen  values  and 
eigen  functions  of  very  large  matrices  (order  of  matrix  =  number  of  vibrational  levels),  so 
that  approximations  are  eventually  introduced. 

Montroll  and  Shuler8  consider  a  "ladder  climbing"  model  in  which  transitions  can  take  place 
only’ between  nearest  neighbor  levels  and  dissociation  only  from  the  highest  level  N 
(dissociation  consists  of  attainment  of  level  N  +  1).  In  the  language  of  the  theory  of  stochastic 
processes,  this  is  the  problem  of  a  one  dimension  random  walk  with  reflecting  (level  zero) 
and  absorbing  (level  N  +  1)  barriers.  The  reciprocal  of  the  "mean  first  passage  time"  for 
the  N'th  level  serves  as  an  average  reaction  rate. 

The  authors  apply  their  general  theory  to  the  case  of  harmonic  oscillators  where  the 

transition  probabilities  are  given  by  the  Landau  Teller  relations  =  [(m  +  1)  bn„l  m  J' 

m  6  1  P  =  p  where  P  is  the  probability  per  collision  of  transition  from  level 

n+l,m  lo  nm  mn 

n  to  level  m,  and  6  is  the  Kroneker  delta.  They  find  that  the  rate  of  chemical  reaction  is 
nm 

not  much  affected  by  deviation  from  vibrational  equilibrium  if  the  activation  energy  e*  is 
larger  than  10  kT. 


A 
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Kim*'*  gives  a  general  mathematical  treatment  of  a  problem  similar  to  that  of  Montroll  and 
Shuler  but  in  which  transitions  are  allowed  between  any  two  levels  (dissociation  only  from 
highest,  level).  He  calculates  the  mean  first  passage  time  for  a  truncated  harmonic 
oscillator  with  an  arbitrary  initial  distribution  of  the  Xj  and  with  a  Boltzmann  initial 
distribution  at  a  vibrational  temperature  ^  T  and  for  a  slightly  anharmordc  oscillator,  where 
nearest  neighbor  and  next-nearest  neighbor  transitions  only  are  allowed. 


Gans*4  presents  the  formal  solution  for  a  set  of  equations  more  general  than  that  of  Kim  in 
that  dissociation  can  take  place  from  any  state  (not  just  the  highest)  and  molecules  may  enter 
the  system  at  any  level. 


Shuler15,  treating  the  same  model  for  dissociation  of  a  truncated  harmonic  oscillator  as 
Montroll  and  Shuler,  uses  their  general  solution  for  the  vibrational  level  distribution  Xj  (t)  *n 


terms  of  Gottlieb  polynomials.  He  specializes  this  to  the  cases  of  a  Boltzmann  and  a  delta 
function  (Xj  (0)  =  6^)  initial  distribution.  Some  numerical  work  for  an  initial  Boltzmann 
distribution  presented  as  a  plot  of  In  Xn  versus  n  shows  very  little  perturbation  of  this 
distribution  for  e  *  AT  =  11  (the  largest  value  considered),  the  population  of  only  the  highest 


level  (n  =  10)  being  appreciably  reduced.  He  also  gives  a  figure  shi/wing  the  effect  of 
dissociation  on  the  rola  ;ation  of  a  system  of  oscillators  all  of  which  are  initially  in  the 
level  n  =  3,  where  the  highest  discrete  level  is  n  =  4.  In  this  work  the  probability  of 
dissociation,  transition  from  level  N  to  "level"  N  +  1,  is  taken  as  (N  +  1)  P^#  i.  e. ,  as 
though  N  +  1  were  the  next  vibrational  level. 


Nikitin16  treats  a  ladder  model  for  dissociation  of  a  harmonic  oscillator  similar  to  the  above 
except  that  probability  of  dissociation  is  taken  as  a  parameter.  Necessary  eigen  values 
and  eigen  vectors  are  found  by  a  perturbation  method,  and  simple  approximate  expressions 
are  derived  for  the  reaction  rate  and  the  distribution  Xr  (t)  that  hold  after  a  certain  induction 
time.  They  find  that  for  small  P  /P,  _  the  distribution  Xn  (t)  differs  little  from  the  Boltzmann 
but  for  large  P  /P1Q  the  equilibrium  distribution  is  disturbed,  so  that  the  probability  of  finding 


a  molecule  in  an  upper  level  is  close  to  zero.  If,  further,  P  /P  »  ?  */kT,  the  reaction 

...  O0  10 

rate  is  independent  of  P  . 

oo 
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In  a  later  paper  Nikitin  and  Sokolov  apply  the  above  model  and  methods  to  an  arbitrary 
oscillator.  They  obtain  for  the  dissociation  rate,  after  an  induction  time,  an  approximate 
equation  in  which  the  transition  probabilities  do  not  appear  and  the  energy  levels  ente  * 
through  Q,  the  vibrational  partition  function.  Before  confronting  this  dissociation  rate  with 
experimental  results,  the  authors  correct  It  for  the  influence  of  rotation  on  the  effective 
dissociation  energy  and  the  Influence  of  vibrational  strte  on  collision  diameter  and  multiply 
by  the  electronic  statistical  weight.  The  equation  can  be  extended  to  polyatomic  molecules, 
and  a  table  of  experimental  and  calculated  values  of  the  ratio  of  pre -exponential  factor  to 
collision  number  shows  excellent  agreement  for  both  diatomic  and  polyatomic  molecules. 

Nikitin  and  Sokolov  point  out  that  for  an  anharmonic  oscillator  dissociation  leads  to  an 
essential  perturbation  of  the  Boltzmann  distribution;  because  the  energy  levels  become 
very  close  together  near  the  dissociation  limit  (see  Figure  l18),  there  are  many  levels  In 
the  energy  range  kT  and  hence  many  levels  with  depleted  population.  This  conclusion  is 

opposite  to  that  which  can  be  drawn  from  Shuler’s  paper  for  a  truncated  harmonic  oscillator 
when  e*/kT  »  1. 

j 9 

Eyring  and  collaborators  give  a  formal  solution  to  a  set  of  transport  equations  which  is 
more  general  than  that  used  by  Kim  in  that  molecules  can  be  injected  at  any  level  and  less 
general  than  those  of  Cans  in  that  only  the  N  +  1  level  is  absorbing.  An  explicit  solution  is 
then  given  for  the  case  of  equally  spaced  levels  where  only  nearest  neighbor  transitions  are 
allowed  and  the  transition  probability  Pfc+1^  is  independent  of  the  level  K  (for  the  harmonic 
oscillator  P^  fc  Is  proportional  to  K  +  1).'  If  all  the  molecules  are  initially  in  the  lowest 
state,  and  no  new  molecules  are  Injected  into  the  system,  it  is  found  that  when  e  */kT»  1 
the  ratio  of  reaction  rate  k  (reciprocal  of  mean  first  passage  time)  to  the  equilibrium  rate 

ke  is  gi  /en  by  k/\  =  1  "  6  T*  where  6  is  the  seeing  between  energy  levels.  (It  was 
also  assumed  that  6/2kT  >  3/N  +  1  or  e  */kT  >  6).  This  identical  result  was  obtained  by 
Montrose  and  Shuler  for  the  truncated  harmonic  oscillator  when  N  »  l. 


> 


Again,  let  all  moleculeo  be  initially  in  the  lowest  level  and  let  molecules  be  injected  into 
this  level  at  a  constant  rata.  After  sufficient  time,  a  steady  state  will  be  set  up,  and  the 
authors  find  that toe  ratio  of  steady  state  to  the  equilibrium  reaction  rate  is  also  given  by 


k  /k  =  1-e 
ss  e 


2. 1.  3  SOME  TREATMENTS  OF  NONEQUI LIBRIUM  REACTIVE  FLOW 

lie  methods  for  treatment  ot  reactive  flews  in  several  papers  are  described  below  in  the 

order  of  the  completeness  with  which  the  authors  consider  various  nonequilibrium  processes. 

2. 1.  3. 1  Finite  Reaction  Rates  Only 

t,  ,20 

c  assumes  that  both  rotation  and  vibration  are  in  equilibrium  with  translation  and  that 
the  approach  to  local  thermodynamic  equilibrium  is  governed  by  the  dissociation  and  re¬ 
combination  rates  in  the  reactions  0+M*0  +  0+  MandN  +M*N  +  N+M. 

Dissociation  rates  kd  are  calculated  by  Logan's  modification  of  simple  collision  theory21, 
in  which  molecules  in  different  vibrational  levels  are  considered  separate  species  and  the 
prol  .bllitv  of  dissociation  by  a  sufficiently  energetic  collision  is  the  same  for  all  levels. 
Recombination  rates  kr  are  calculated  from  kj/k.  -  K  where  Ke  is  the  equilibrium  constant. 
Die  use  of  this  equation  is  consistent  with  the  assumption  of  vibrational  equilibrium,  but  its 
use  is  not  justified  where  vibrational  equilibrium  does  not  hold*7,22*23. 


2*1*3*2  gjnite  Reaction  Rate  Plus  Vibrational  Relaxation 

An  earlier  paper  of  Lick's24  is  very  similar  to  the  above.  Here,  however,  vibration  is  not 
assumed  to  be  in  equilibrium  with  translation,  and  the  vibrational  energy  is  calculated  from 
the  Landau  Teller  equation  for  the  relaxation  of  a  set  of  harmonic  oscillators:  dE  /dt  = 
fv  "  VT  where  \  is  the  equilibrium  vibrational  energy  at  the  temperature  T. 

2. 1.  3.  3  Interaction  of  Dissociation  and  Vibration  (First  Model} 

Duff  and  Davidson26  consider  the  coupling  between  vibrational  excitation  and  chemical 

reaction  by  adopting  the  following  model:  The  C>2  molecules  (or  N  )  are  considered  to 
consist  of  two  classes: 


02  -  those  In  the  ground  vibrational  state 


* 

Og  -  all  those  in  excited  vibrational  states 


Dissociation  was  assumed  to  take  place  only  from  excited  states  and  the  following  reactions 
assumed  to  occur: 

o2  +  m^o2  +  m  n+o2*no  +  o 


nJ  +  M  ss  N*  +  M 
2  2 


$ 

N+O  si  NO  +  O 

Li 


HfL 

O+M^O+o  +  m 

Li 


3|t 

N+MssN+N  +  M 

Li 


NO  +  M2N+0  +  M 

Numerous  kinetic  assumptions  had  to  be  made  in  the  use  of  these  equations. 

*  * 

The  02  (and  Ng)  is  assumed  to  be  in  equilibrium  with  the  translational  heat  bath;  i.  e. ,  the 
relative  population  of  the  excited  vibrational  levels  are  given  by  a  Boltzmann  distribution  at 
the  translational  temperature.  This  assumption  was  made  for  convenience  in  calculation;  it 
avoids  the  introduction  of  a  vibrational  temperature. 


The  necessary  thermodynamic  functions  for  the  two  classes  of  molecules  are  obtained  by 
splitting  the  partition  function  of  the  harmonic  oscillator  into  a  part  corresponding  to  the 
ground  state  and  a  part  corresponding  to  all  the  excited  states. 


igteraction  oi  Dissociation  and  Vibration  (Second  Model) 

A  more  sophisticated  model  for  the  Interaction  of  dissociation  and  vibration  has  been 
introduced  by  Hammerling,  Teare,  and  Kivel22  and  subsequently  used  by  other  authors 
Their  method  is  based  on  the  assumption  that  dissociations  can  take  place  from  any 
vibrational  level  and  that  t'  e  efficiency  of  collision  in  producing  dissociation  from  a  giv 
level  is  the  same  for  all  levels.  The  dissociation  rate  is  then  taken  to  be  proportional  1 
the  sum  over  all  levels  of  the  product  of  the  probability  of  being  a  level  and  the  fraction 
molecules  having  translational  energy  D  -E  or  greater: 


where 


fraction  of  molecules  in  n’th  level 


total  number  of  vibrational  states 


dissociation 


local  heat  bath  temperature 


To  use  this  formula  requires  knowledge  of  the  X  under  conditions  of  vibrational  non- 
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equilibrium.  For  this  purpose  the  authors  make  use  of  the  results  of  Montroll  and  Shuler  , 

who  showed  that  for  a  dilute  solution  of  harmonic  oscillators  in  a  constant  temperature  heat 

bath,  an  initial  Boltzmann  distribution  will  relax  to  a  fins  1  Boltzmann  distribution  through  a 

continuous  sequence  of  Boltzmann  distributions  characterized  by  a  parameter  T,,»  the  so 
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called  "vibrational  temperature.  "  In  fact,  Bazley,  Montroll,  and  Rubin  show  this  to  be 
nearly  true  also  for  anharmonic  oscillators  with  a  Morse  potential.  They  present  the 
following  equation  for  the  ratio  of  instantaneous  to  equilibrium  dissociation  rate: 


where 


vibration  constant 


©v  can  be  obtained  from  the  equation  for  the  vibrational  energy  relaxation  of  harmonic 
oscillators: 


where 


average  vibrational  energy  per  modecule 
vibrational  relaxation  time 


and  the  relations 


(e)eq  -  hcwe/(e°  -  1) 

0 


=  hccu  /  (e  -  1) 


(2. 1.3-4) 
(2. 1.3-5) 


2. 1.3. 5  Incorporation  of  Effect  of  Dissociation  on  Vibrational  Energy 

Dissociation  does  not  take  place  at  an  equal  rate  from  all  levels;  highly  excited  molecules 

having  vibrational  energies  close  to  are  most  easily  dissociated.  It  is  clear,  therefore, 

that  the  dissociation  process  reduces  the  average  energy  remaining  in  the  molecules. 
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Treanor  and  Marrone  adopt  the  method  of  Hammerling,  Teare,  and  Kivel  but  take  this 
effect  into  account.  They  assume  that  dissociation  does  not  appreciably  perturb  the 
Boltzmann  distribution  of  the  vibrational  levels  but  modify  the  equation  for  relaxation  of 
vibrational  energy  to  include  loss  of  energy  by  dissociation  and  gain  by  recombination: 


.  (e)  -  e 

de  _  eq 

dt  r 


[I  (T,  Ty)  -  e] 

[X] 


M?) (HP) 


where 


[x]  =  number  density  of  molecules 

f,  r  =  forward  (dissociation),  rearward  (recombination) 

E(T,  T  )  =  average  energy  lost  from  vibration  in  a  single  dissociation 
G(T)  =  average  energy  gained  by  vibration  in  a  single  recombination 


From  the  assumptions  concerning  dissociation  made  above  (Section  2. 1. 3. 4),  it  can  be  shown 
that  for  a  truncated  harmonic  oscillator 


At  equilibrium  de  /dt-*0,  T+  T,  e  -*  (e)^,  and  (d  [ X  ] / dt)r*  (d  [X]/dt)^  it  follows  that 

G(T)  =  T  E(T,  Ty)  (2. 1.3-8) 

V 

=  \  (N  -  1)  w 

Calculations  made  by  the  authors,  using  the  method  of  Hammer  ling,  Teare,  and  Kivel  without 
modification,  of  the  temperature  distribution  behind  a  normal  shock  in  oxygen  or  in  nitrogen 
show  that  the  vibrational  temperature  overshoots  the  translational  temperature  T  at  the 
higher  flow  speeds  (a  similar  result  for  O  was  obtained  in  Reference  4).  Incorporation  of 
the  modified  vibrational  relaxation  equation  results  in  higher  translational  and  lower 
vibrational  temperatures,  completely  eliminating  the  vibrational  temperature  overshoot. 


2.1.4  DISCUSSION 

As  mentioned  above,  Montroll  and  Shuler's  work  on  the  relaxation  of  a  system  of  harmonic 
27 

oscillators  ,  upon  which  Hammerling  et.  al.  based  the  assumption  of  a  Maxwell-Boltzmann 

distribution  of  vibrational  levels  at  all  times,  assumes  that  the  oscillators  are  contained  in  a 

large  excess  of  inert  gas  that  acts  as  a  constant  temperature  heat  bath.  Although  this 

condition  holds  for  many  shock  tube  studies,  where  the  gas  under  investigation  is  mixed  with 

a  much  larger  amount  of  argon,  it  will  not  hold  for  the  flow  of  air  about  a  reentry  body. 

Even  here,  however,  if  the  temperature  is  not  too  high,  the  large  excess  of  unexcited 

oscillators  may  be  considered  as  serving  as  the  heat  bath  for  the  small  fraction  of  excited 
27 

oscillators  .  Of  course,  the  heat  bath  temperature  is  far  from  constant  in  hypersonic  flows. 

The  essential  point,  however,  is  that  Reference  27  treats  relaxation  only;  dissociation  was 
not  included.  And  it  can  be  concluded  from  the  material  of  Section  2. 1. 2  that  for  anharmonic 
oscillators  (real  molecules)  dissociation  may  deplete  the  population  of  upper  levels  and  thus 
considerably  perturb  the  Maxwell-Boltzmann  distribution. 
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The  eqi  ition  used  for  dissociation  rate  is  based  on  the  assumption  that  the  efficiency  of 

sufficiently  energetic  collisions  in  producing  dissociation  from  a  given  vibrational  level  is 

the  same  for  all  levels.  This  assumption  was  made  in  the  absence  of  specific  knowledge  of 

these  efficiencies,  but  actually  it  is  not  consistent  with  the  well-known  difficulty  of  converting 

30 

translational  energy  into  vibrational  energy  .  That  the  probability  of  transfer  of  trans¬ 
lational  energy  varies  inversely  as  the  energy  to  be  transferred  is  shown  by  the  anharmonic 

31  23 

(Morse)  oscillator  calculations  of  Heaps  and  Herzberg  ,  Pritchard  ,  and  Bauer  and 
32 

Cummings  . 

The  situation  for  transitions  between  bound  states  should  also  hold  for  transitions  from 

bound  states  to  the  continuum  (i.  e. ,  dissociation)  and  suggests  that  the  efficiency  of  collisions 

in  producing  dissociation  and  the  rate  of  dissociation  from  any  but  the  highest  levels  will  be 

insignificant.  There  is  very  little  available  in  the  way  of  calculations  against  which  to  check 

23 

this  conclusion;  however,  some  crude  calculations  based  on  Pritchard's  work  show  that 
thr  three  highest  states  make  about  equal  contributions  to  the  total  dissociation  rate,  the 
contributions  of  others  being  negligible.  In  calculating  matrix  elements  for  the  collisional 
perturbation  potential,  Pritchard  uses  the  wave  fu  action  for  the  highest  bound  state  (M  =  14) 
in  lieu  of  the  wave  function  for  the  continuum. 

A  much  more  complete  and  general  analysis  of  dissociation  rates  from  various  vibration 
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levels  has  been  made  by  Herzfeld  .  He  proceeds  from  his  work  c  a  relaxation  and  extend^ 

34  35 

some  results  of  Lennard-Jones  and  Strachan  and  Strachan  in  order  to  calculate  matrix 
elements  involving  wave  functions  for  the  continuum.  Herzfeld's  final  result  is  only  formal 
in  that  there  still  remains  to  be  carried  out  the  integration  of  a  series  (finite)  of  complicated 
terms  involving  gamma  functions  of  complex  arguments;  but  he  shows  how  approximate 
calculations  can  be  made  for  dissociation  from  the  ground  state  and  from  states  near  the  top 
of  the  Morse  curve.  Numerical  work  on  O  in  Herzfeld's  paper  shows  that  only  the  upper- 
most  four  or  five  levels  contribute  to  the  dissociation  rate. 
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Although,  as  is  evident  from  the  above,  the  assumptions  underlying  the  model  of  Hammerling 
et.  al.  may  deviate  considerably  from  reality,  the  model  was  forced  to  give  correct  results 
at  equilibrium  and  hence  reasonable  approximations  near  equilibrium.  This  was  accomplished 
by  using  it  for  calculating  not  absolute  reaction  rate  but  the  ratio  of  dissociation  rate  under 
(vibrational)  nonequilibrium  to  the  equilibrium  rate. 

It  would  be  possible  to  construct  a  somewhat  more  realistic  model  in  which  dissociation  takes 
place  only  from  a  few  highest  levels,  perhaps  with  tne  additional  assumption  of  equal  prob¬ 
ability  of  dissociation  from  any  of  these  levels  under  sufficiently  energetic  collisions.  This 
would  require  calculation  of  the  time  history  of  the  population  of  these  levels  as  well  as  the 
time  history  of  the  total  vibrational  energy.  These  quantities  could  be  obtained  by  setting  up 
a  ladder  model  incorporating  dissociation  from  several  upper  states,  calculating  transition 
and  dissociation  rates  from  Herzfeld's  work,  and  solving  the  resulting  very  numerous 
transport  equations  by  the  matrix  methods  of  the  references  given  in  Section  2.1.3. 

On  the  conditk  .hat  the  coupled  vibrational-dissociational  problem  need  be  considered  only 
for  the  ground  electronic  state  (which  is  certainly  questionable  for  hypersonic  shock  layer 
applications),  this  type  of  solution  should  be  possible  for  simple  systems  such  as  oxygen, 
which  has  between  21  and  53  vibrational  states  in  the  ground  electronic  state,  or  nitrogen, 
which  has  between  30  and  70  vibrational  states  in  the  ground  electronic  state.  For  air, 
however,  with  the  addition  of  between  25  and  60  vibrational  states  for  nitric  oxide  (which 
comes  to  a  total  of  between  76  and  183  vibrational  states  that  must  be  accounted  for)  this 
type  of  solution  becomes  almost  impossible.  Considering  this  plus  the  fact  that  as  yet  there 
is  no  method  available  to  account  for  vibrational  interaction  in  reactions  such  as: 

NO  +  O  tt  02  +  N 

N  +  O  tt  NO  +  N 


2NO  **  N  +  0„ 

2  2 


leads  to  the  conclusion  that  this  type  of  solution  fcr  the  air  system  is,  indeed,  at  present 
impossible. 
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Subsequent  to  Treanor  and  Marrone’s  work  ,  Marrone  published  the  development  of  a 
preferential  model  for  vibration  dissociation  coupling  which  provides  for  variable  probability 
of  dissociation  from  different  vibrational  levels  in  the  form 

F  (v)  =  exp  -  [D  -  e  (v)  ]/kV  (2. 1.4-1) 


where  v  denote  the  vibrational  level;  D,  the  dissociation  energy  from  the  ground  (v  =  O) 
level;  e  (v)f  the  energy  of  the  vibrational  level;  k,  Boltzmann’s  constant;  and  V,  an 
arbitrary  parameter.  The  energy  levels  are  obtained  from 


2  3 

e  (v)  =  (v  +  1/2)  -  we  Xe  tv  +  1/2)  +  (v  +  1/2) 


(2. 1.4-2) 


in  which  the  coefficients  are  determined  from  spectroscopic  data.  A  Boltzmann  distribution 
is  used  for  the  vibrational  level  populations.  Although  certainly  limited  by  this  last  condition 
as  by  the  use  of  an  arbitrary  parameter  to  prescribe  transition  probabilities,  Marrone’s 
preferential  model  for  coupled  vibration-dissociation  is  the  best  available  at  present. 

2. 1. 5  SUPERPOSmONAL  METHOD  FOR  INVESTIGATING  VIBRATIONAL  NONEQUILIBRIUM 

OQ 

Landau  and  Teller*^ J  have  shown  that  the  Master  Equation  for  a  system  of  harmonic 
oscillators  reduces  to  the  familiar  vibrational  relaxation  equation: 


(2. 1.5-1) 
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in  which  R  is  the  universal  gas  constant,  8^  is  the  characteristic  vibrational  temperature; 

T,  the  kinetic  or  gas  temperature;  and  T^,  the  vibrational  temperature.  For  applications  in 
which  the  relaxation  time  and  the  equilibrium  vibrational  energy  are  implicit  functions  of 
time,  the  general  solution  of  equation  1  is: 


dt*  +  const 


where  t*  and  t**  are  dummy  variables  for  time.  Applying  the  boundary  condition  at  t 
€  *  €  allows  the  integrals  to  be  put  into  the  definite  form: 


For  steady  flow  the  distance  along  a  streamline,  s,  and  the  ordered  velocity,  V,  are  related 


so  that  the  integrals  in  Equation  2.  i.  5-6  may  be  transformed  into: 


‘f-4  -^vr22  exp(/  ^r)  H  exp(-/  ^)+%„  p-1-5-7) 


Using  Equations  2. 1. 5-2  and  2. 1.  5-3  the  vibrational  energy  may  be  rewritten  as: 


H,v  !•(  VT  [  exp  (Bv/T)  -  l  -  exp  te^/T)  -  !  ]  e*p(f0  ^r)jexp(-f0 


ds*  j 

Vr/ 


exP  <VTvo>  - 1 


(2. 1.5-8) 


and  the  local  vibrational  temperature  along  the  streamline  is  given  by 


T  =  9  /In  (1  +  R  0  /fe  ) 

V  v  V  V 


(2. 1.  5-9) 


By  superimposing  the  solution  of  Equations  2. 1. 5-8  and  2. 1. 5-9  over  the  thermal  equilibrium 
calculations  presented  in  the  first  quarter  report,  the  relative  importance  of  the  vibrational 
nonequilibrium  problem  can  be  assessed. 

2.2  FUTURE  FLOW  FIELD  WORK 

The  literature  survey  will  be  completed  and  presented  in  the  third  quarter  report. 

f  ..v 

The  results  of  the  superimposed  vibrational  nonequilibrium  calculations  obtained  by  the 
method  described  above  and  a  critique  on  the  relative  Importance  of  this  problem  will  be 

I n  flio  _ a 
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Shock  layer  stagnation  line  profiles  from  the  coupled  influscid-boundary  layer  chemical 
nonequilibrium  solution  as  well  as  any  nonequilibrium  boundary  layer  results  which  may  be 
available  will  be  reported  at  the  end  of  the  third  quarter. 
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SECTION  3 

REENTRY  EFFECTS  ON  ANTENNA  BREAKDOWN 


The  subject  of  this  section  is  the  effects  of  reentry  on  voltage  breakdown  of  antennas,  with 
emphasis  on  the  reentry  environment  associated  with  slender  cone  bodies  at  hypersonic 
speeds.  The  discussion  will  depend  entirely  on  ideas  and  data  taken  from  written  reports  and 
verbal  communications  with  a  number  of  researchers  in  the  field,  in  the  manner  of  a  survey 
of  present  knowledge. 


For  conical  vehicle  re-entry,  three  distinct  flow  field  regions  may  be  important  in  antenna 
breakdown:  the  ambient  free  stream,  the  inviscid  shock  layer,  and  the  boundary  layer.  With 
the  exception  of  relatively  high  altitudes  where  these  regions  are  not  necessarily  distinct, 
the  reentry  environmental  effects  on  breakdown  will  be  different  for  each  of  the  three  regions. 
In  fact,  most  workers  assume  that  a  convenient  separation  can  be  made  such  that  the  effects 
of  each  region  can  be  treated  independently.  At  any  rate,  the  different  basic  processes  which 
control  breakdown  can  be  identified  for  each  of  the  three  flow  field  regions,  for  example,  as 
done  by  Reilly  and  numerous  other  authors. 

The  free  stream  ahead  of  the  shock  layer  is  characterized  essentially  only  by  the  fact  that 
the  air  is  moving  relative  to  the  transmitting  antenna.  The  effects  of  altitude  on  the  state  or 
composition  of  the  ambient  air  probably  are  not  significant  in  breakdown;  and  even  if  they 
were,  it  would  be  a  relatively  simple  task  to  account  for  them.  Hence,  the  only  important 
effect  in  the  free  stream  is  convection. 


The  inviscid  flow  region  is  characterized  by  a  stfite  of  elevated  temperature  and  increased 


particle  number  density  relative  to  the  free  stream.  On  a  slender  conical  vehicle,  ionization 
is  probably  insignificant  in  this  region,  but  the  detailed  thermodynamic  state  and  chemical 
composition  of  the  shock  heated  air  maj  be  quite  important.  At  any  rate,  the  effects  of  the 
inviscid  flow  region  are  felt  in  terms  of  convection,  increased  air  density,  and  elevated 


e  oun  ary  layer  is  characterized  by  eyen  higher  temperatures  produced  by  viscous 
dissipation.  The  presence  of  significant  levels  of  thermal  ionization  is  thus  probable. 
Particle  number  density  and  convection  are  not  now  so  easy  to  generalize  because  of  tl 
strong  gradients  of  these  quantities  which  exist  here.  In  addition,  the  presence  of  chei 
contamination  of  the  air  by  products  of  ablation  will  be  important.  Thus,  convection,  > 
density,  temperature,  ionization,  and  chemical  additives,  all  in  the  context  of  strong 
gradients,  must  be  considered  in  the  boundarv  laver 


where  D  is  the  diffusion  coefficient.  Under  the  aesuMptlon  of  lncon>Pre..1bto  Eow,  tte 
continuity  equation  becomes 


where  v  is  the  net  ionization 
has  been  obtained  by  Chown46 


frequency.  For  the  case  of  constant  D,  an  approximate  solution 
using  an  effective  diffusion  length  A. 


3.2  GAS  DENSITY  EFFECTS 


In  voltage  breakdown  experiments  under  standard  laboratory  conditions,  the  practice  is  to 
express  the  state  of  the  gas  being  tested  in  terms  of  pressure,  since  this  is  a  directly 
measurable  quantity.  However,  pressures  in  the  reentry  induced  shock  layer  are  often 
orders  of  magnitude  higher  than  in  the  free  stream.  Also,  because  of  increased  gas  tempera¬ 
ture,  shock  layer  density  is  not  generally  the  same  as  in  the  free  stream.  MacDonald47  has 
noted  that  the  physically  important  parameter  for  breakdown  is  ihe  gas  density  rather  than 
pressure,  since  it  is  the  number  of  gas  particles  per  unit  volume  which  determines  the  rates 
for  all  the  processes  which  control  breakdown.  Thus  the  use  of  either  the  free  stream 
pressure  or  the  shock  layer  pressure  in  predicting  breakdown  at  a  given  altitude  would  give 
erroneous  results.  The  distribution  of  gas  density  in  the  high  field  intensity  region  is  the 
correct  determining  factor. 


In  order  to  give  a  quantity  which  is  commensurable  with  the  classically  used  parameter  of 
pressure,  most  authors  define  a  reduced  pressure  p*  Reilly39  defines  the  reduced  pressui 
in  terms  of  the  gas  number  density  N  in  particles  per  cubic  centimeter. 


Light  and  Taylor  define  the  reduced  pressure  in  an  equivalent  way  in  terms  of  temperature 
T  in  °K: 


wnerc  p  is  the  local  pressure.  These  correctirus  have  the  effect  of  shifting  the  curve  of 
breakdown  field  versus  altitude  about  30, 000  feot  higher  in  altitude  at  velocities  of  the  order 
of  22, 000  ft/sec,  according  to  MacDonald47. 


These  effects  have  been  observed  in  reentry  flight  experiments.  Bisbing49  showed  that  the 
observed  period  of  breakdown  of  a  slot  antenna  on  a  conical  vehicle  agreed  with  laboratory 
data  in  cold  air  when  the  air  density  corresponding  to  the  inviscid  flow  region  was  used  as  the 
independent  variable.  Also,  Nanevicz50  showed  that  when  the  VHF  breakdown  data  from 
Nike-^ajun  flights  is  corrected  for  local  density,  it  agrees  roughly  with  laboratory  measure¬ 
ments  in  cold  air. 

3.3  HIGH  TEMPERATURE  EFFECTS 

The  effects  of  elevated  gas  temperature  on  voltage  breakdown  have  been  investigated  in  terms 
oi  dc  breakdown  at  a  relatively  early  date**4.  It  was  found  that  up  to  temperatures  of  1400°K 
no  deviation  from  cold  gas  results  could  be  observed52.  The  investigation  of  this  phenomenon 
by  the  use  of  a  shock  tube  in  order  to  produce  higher  temperatures  was  carried  out  by 
Sharbaugh  et.  al.  for  nitrogen  .  These  investigators  found  that  significant  departures  from 
cold  gas  breakdown  occur  at  temperatures  of  2000°K  and  that  at  4000°K  the  breakdown  potential 
is  half  the  cold  gas  value.  Their  analysis  of  this  effect  was  based  on  a  theory54  of  space- 
charge  alteration  of  the  field  by  ionization  in  the  shock  tube,  although  certain  anomalies  were 
present,  which  they  ascribed  to  ionization  relaxation  in  the  shock  tube.  The  same  investi¬ 
gators  carried  out  further  experiments  in  both  nitrogen  and  air,  in  which  the  electron  density 
was  also  measured  .  These  data  showed  an  apparent  correlation  between  the  reduction  in 
breakdown  voltage  and  the  electron  density. 

Notwithstanding  the  fact  that  dc  breakdown  is  affected  by  surface  phenomena  and  space- 
charge  effects,  it  may  be  that  the  above  results  provide  an  indication  of  effects  of  tempera¬ 
ture  on  basic  processes  of  gaseous  breakdown.  In  fact,  it  is  apparent  that  they  are  in 
qualitative  agreement  with  more  recent  work  on  effects  of  high  gas  temperature  on  microwave 
breakdown5  ’ 57 » 48.  Light  and  Taylor48  have  attempted  to  explain  the  lowering  of  the  break¬ 
down  field  at  high  temperatures  in  terms  of  the  effect  of  vibrationally  excited  states  of  oxygen. 

In  their  analysis  of  the  data,  they  postulated  that  the  curve  of  ionization  frequency  versus 
field  strength  is  shifted  to  lower  field  strength  values  by  the  mechanisms  associated  with 
electron  heating.  The  mechanism  which  they  use  to  explain  the  data  is  inelastic  collisions 


between  electrons  and  oxygen  molecules.  When  the  oxygen  is  vibrationally  excited  in  high 
temperature  air,  the  electrons  are  more  easily  heated  by  the  electric  field  due  to  the  lack  of 
important  electron  energy  transfer  to  ground  state  molecules.  A  formula  for  this  effect  is 
postulated  as  follows. 


vi  _  f  /  v2  ] 

P*  V  (T,  p*)  \  p*  J  j  (3. 3-i) 


where  Vj  is  the  net  ionization  frequency,  p*  is  reduced  pressure,  T  is  temperature,  T  is 
273  K,  and  Eg  is  effective  field.  The  effect  of  oxygen  excitation  is  placed  in  the  form 


ns l  p*)  _ 

V(T0) 


kT\  XJTj_  p*) 

/  X<V 


(3. 3-2) 


where  A  is  the  minimum  energy  required  to  excite  a  vibrational  state  (0. 2  ev  for  oxygen),  k 
is  Boltzmann's  constant,  and  X  is  the  mole  fraction  of  oxygen  molecules. 

The  equations  of  Light  and  Taylor  predict  that  the  amount  of  molecular  oxygen  present  would 

have  a  strong  effect  on  the  ionization  frequency  as  a  function  of  gas  temperature.  However 
57 

Thompson  ,  in  an  arc  channel  experiment,  observed  no  apparent  difference  between  pure 
nitrogen  and  mixtures  of  nitrogen  and  oxygen.  This  result  would  tend  to  place  the  hypothesis 
of  Light  and  Taylor  in  doubt,  except  that  the  arc  channel  experiment  is  probably  more  sensitive 
to  influences  from  other  phenomena.  At  any  rate,  the  hypothesis  that  oxygen  excitation  alone 
is  import"  t  is  an  obvic;  s  oversimplification  in  that  the  absence  of  oxygen  would  certainly 
require  some  other  mechanism  of  electron  energy  transfer  to  dominate.  Note  also  that  in  the 
equation  of  Light  and  Taylor  oxygen  dissociation  has  a  much  greater  effect  in  terms  of  tempera¬ 
ture  than  does  excitation,  since  dissociation  itself  is  a  strong  function  of  temperature. 


Aside  from  the  effects  of  ionization,  which  are  discussed  elsewhere,  various  authors  have 
identified  other  possibly  important  effects  of  high  temperature  on  breakdown,  Reilly39 
suggests  that  molecular  excitation  promotes  ionization  by  effectively  lowering  the  energy 
required  to  ionize  the  molecule.  Epstein58  suggests  that  the  gas  composition  is  important 
in  terms  of  the  presence  of  nitric  oxide,  which  ionizes  more  readily  than  cold  air.  Epstein 
also  notes  that  attachment  and  detachment  rates  increase  with  increasing  air  temperature. 


3.4  IONIZATION  EFFECTS 

A  significant  amount  of  thermal  ionization  is  known  to  exist  in  the  shock-heated  air  surround¬ 
ing  a  reentry  vehicle  .  The  resulting  high  initial  electron  density  nQ  which  exists  before 
the  electric  field  is  turned  on  may  influence  the  rates  of  several  of  the  electron  production/ 
loss  mechanisms.  The  following  discussion  explores  three  possibxe  changes  in  the  breakdown 
process  due  to  high  initial  electron  density  as  considered  in  the  literature. 


The  first  and  probably  most  important  effect  of  large  electron  densities  is  to  change  diffusion 
from  a  free  to  an  ambipolar  process.  That  is,  with  the  buildup  of  a  sizable  election  space 
charge,  the  electron  and  molecular  diffusion  processes  become  inextricably  linked.  The 
result  is  to  slow  down  the  diffusion  of  electrons  to  the  molecular  rate  and  thereby  correspond¬ 
ingly  reduce  diffusion  as  an  electron  loss  mechanism  (lowering  the  breakdown  field).  Allis 
and  Rose  have  indicated  that  the  transition  from  free  to  ambipolar  diffusion  begins  where 
the  Debye  length  becomes  comparable  to  the  diffusion  length.  If  we  make  the  approximation 
that  the  diffusion  length  is  approximately  half  the  applied  wavelength,  then  we  must  consider 
ambipolar  diffusion  for  electron  densities  greater  than  about  106  electrons/cm3  for  signal 
frequencies  on  the  order  of  100  megahertz.  Such  electron  densities  clearly  occur  throughout 
reentry  in  at  least  some  part  of  the  flow  field59  and  therefore  require  some  modification  of 
the  diffusion  coefficient.  Whitmer  and  MacDonald01  have  indicated  that  the  free  diffusion 
coefficient  D  is  related  to  the  ambipolar  diffusion  coefficient  D?  by  the  equation  D  =2DM  /M 
where  p+  and  \x_  are  the  mobilities  of  ions  and  electrons  respectively.  They  show  from 
measurements  of  mobilities  that  the  ratio  D/D&  is  approximately  40  in  air.  Now,  using 
MacDonald's  *  expression  for  free  diffusion  in  terms  of  the  effective  applied  field, 


ill 
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DP  -  (29  ♦  0. 8  E  /P)  104  torr  cm2/seo,  approximate  value,  for  purely  free  or  purely  ambl- 
polar  diffusion  ma^  be  had.  Chown  et  ah  40  suggest  that  the  highly  nonuniform  electron 
density  In  the  flow  field  of  a  reentering  vehicle  makes  it  extremely  difficult  to  determine 
whether  diffusion  will  b»  free  or  ambipolar  or  essentially  a  combination  of  the  two.  In 
addition,  the  fact  that  the  transition  from  free  to  ambipolar  diffusion  may  take  place  over 
several  orders  of  magnitude  change  in  electron  density  makes  it  desirable  to  define  an 
effective  diffusion  coefficient  D,  which  will  take  into  account  both  processes.  Bachynsld 
et.  aL  62  have  defined  such  an  effective  diffusion  coefficient  =e  =  Da  Cl  +  M_  (N-n)/(p+N  +  p_n)] 
(where  N  and  n  are  the  ion  and  electron  number  densities  respectively)  which  can  easily  be 

established  from  the  flow  field  properties. 

The  second  effect  of  high  Initial  electron  densities  Nq  is  to  shorten  the  time  required  to  reach 
critical  electron  density  (n  - 1.24  x  lO’V,  where  f  is  the  signal  frequency)  and  considerably 
reduce  the  breakdown  field  strength.  For  vbf  signals  on  the  order  of  100  megahertz,  the 
critical  density  Is  about  108  electrons/cm3.  which  may  easily  occur  purely  from  thermal 
generation  during  parts  of  a  hypersonic  flight59.  Chown  et.  al.  state  thal  although  the 
details  of  propagation  In  such  an  overdense  plasma  are  net  well  understood,  there  is  experi¬ 
mental  evidence  to  suggest  t  rat  once  the  electric  field  Is  turned  on  the  plasma  will  grow  u> 
thickness  such  that  the  fields  penetrating  beyond  it  are  just  below  the  levels  required  to 
-w.i.  the  discharge  under  the  existing  focal  conditions  of  pressure,  frequency,  etc.  In 
determining  'he  praotibility  of  ECM  system  designs,  then,  It  will  be  Important  to  establish 
over  what  portion  of  the  trajectory  It  can  be  ejected  that  the  transmission  will  be  thermally 
blacked  out  (independent  of  field  strength)  and  just  how  much  the  propagated  signal  will  be 
degraded  by  the  overdense  plasma. 

Finally,  we  would  like  to  consider  the  effect  of  high  electron  densib  on  the  recombination  loss 
rate.  Chown  et.  al. 40  renark  that  there  is  considerable  uncertainty  as  to  the  lower  limit  of 
electron  density  above  which  reeosnbination  is  important,  especially  in  high- temperature  air. 
They  note  that  predictions  have  ranged  anywhere  from  109  to  1012  electrons/cm  and  abo\e. 
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ReMy  considers  the  dissociative  and  the  two  body  olectron-ion  recombination  rates  as  reported 
by  Brown  ,  Massey  and  Burhop  3  over  a  ,*ange  of  electron  densities  from  106  to  1012/cm 
at  NTP.  He  shows  that  the  two-body  dissociative  recombination  rate  remains  several  orders 
of  magnitude  below  the  rate  of  two- body  attachment  for  electron  densities  ranging  up  to 
10  /cm  ,  and  therefore  two- body  recombination  can  be  considered  negligible  at  NTP.  (For 
electron  densities  greater  than  10i2/cm3,  thermal  blackout  will  be  the  predominant  effect.) 

He  likewise  conclude,'  that  since  the  three-body  rate  of  recombination  and  attachment  are  an 
order  of  magnitude  slower  than  two-i»dy  recombination  at  NTP,  they  can  be  safely  ignored. 

For  high-temperature  air,  Beilly  quotes  the  results  of  Lin  and  Tear64,  giving  two-body 
recombination  a  Tf3  2  dependence,  and  the  work  of  Massey  and  Burhop63,  which  shows  a 
:e  dependence  for  the  three-body  rate.  Assuming  these  inverse  temperature  dependencies, 
he  concludes  that  two-body  attachment  dominates  both  two-  and  three-body  recombination  for 
electron  energies  above  1  ev  and  electron  densities  bekrv  l012/cm3  39. 

3.5  EFFECTS  OF  CHEMICAL  ADDITIVES 

The  effects  of  chemtcsl  additives  lo  the  boundary  layer  are  important  because  of  the  natural 

effects  of  ablation  as  well  as  in  terms  of  artificial  breakdown  alleviation  techniques.  In  fact. 

most  of  the  work  done  on  this  subject  appears  to  be  motivated  by  the  search  for  alleviation 
methods. 

Thompson^  has  measured  the  effects  of  CO,  CF^  and  SFg  on  electron  density  and  antenna 
breakdown  in  an  arc  channel.  The  first  two  of  these  materials  are  products  of  ablation  for 
carbon  and  Teflon  heat  shields,  respectively;  and  the  last  is  a  popular  electrophilic  quenchant 
candidate.  Twenty  percent  by  volume  of  each  of  there  materials  was  added  to  the  flow,  and 
experiments  were  conducted  at  temperatures  between  3000°  and  5000°K.  Carbon  monoxide 
was  found  to  have  no  effect  on  either  electron  density  or  breakdown.  CF4  reduced  the  electron 
density  by  a  factor  of  2  to  3  and  increased  the  breakdown  power  by  a  factor  of  3.  SF  .  reduced 
the  electron  density  by  more  than  a  factor  of  3  and  increased  the  breakdown  power  by  a  factor 

of  2.  The  authors  attribute  the  breakdown  alleviation  by  the  fluorides  to  the  additional  electron 
attachment  provided  by  fluorine  and  fluoride  radicals. 


Space  Sciences,  Inc.  of  Waltham,  Massachusetts,  has  been  conducting  analytical  and  ejqjerl- 
mental  studies  on  alleviaticn  of  plasma  sheath  and  antenna  breakdown  effects.  They  have  been 
evaluating  materials  for  the  ability  to  reduce  electron  concentration  and  to  increase  break¬ 
down  power  level  in  thermal  plasmas.  The  primary  alleviation  mechanism  studied  is 
attachment  of  electrons  by  electronegative  gases.  The  plasma  simulation  facility  used  in 
these  studies  is  an  argon  arc  jet.  Conditions  simulated  include  atmospheric  pressure  down 
to  13  torr  ana  temperatures  of  3000  K  to  5000  K.  Diagnostics  are  conducted  by  means  of 
K-band  transmission  through  the  test  section.  Additives  tested  include  SF  ,  BF  ,  N  ,  Freon 
116,  air,  and  argon. 

Space  Sciences  ha3  concluded  that  materials  which  are  good  attachers  in  glow  discharges 
often  do  not  work  in  thermal  plasmas  because  of  dissociation.  A  figure  of  merit  for  attachers 
would  appear  to  depend  on  the  molecular  weight,  the  electron  attachment  cross  section 
(iiioluding  its  energy  dependence),  and  the  resistance  to  tnermal  dissociation.  Highly  elecrr ' 
negative  materials  have  higher  attachment  cross  sections  in  the  molecular  form  than  in  the 
atomic  form.  Survival  lifetimes  of  molecular  additives  in  thermal  plasmas  are  thus  important 
in  determining  their  effectiveness  as  plasma  alleviative  substances. 

The  idea!  additive  material  for  alleviation  of  both  ionization  and  breakdown  is  pictured  by 
Space  Sciences  as  one  made  up  of  "giant  molecules, "  which  would  actually  be  extremely 
small  solid  particles  or  liquid  drops.  These  particles  would  have  maximum  survival  times 
as  well  as  providing  surface  recombination  as  a  mechanism  for  electron  -’emoval.  Refractory 
solids  would  have  to  be  avoided  because  of  thermionic  emission  and  Safca^  angmuir  ionization. 

MIT  Aerophysics  Laboratory  has  been  conducting  experiments  in  a  wind  tunnel  under  the 
direction  of  E.  E.  Covert45.  In  this  experiment  a  supersonic  arc  jet  is  blown  over  the  nose 


of  a  cone  model  in  a  Mach  4  wind  tunnel.  Breakdown  of  an  X-band  slot  on  the  cone  was 
determined  under  conditions  of  injection  of  SFg.  They  have  observed  that  the  electron  density 
was  reduced  by  an  order  of  magnitude,  but  that  the  breakdown  was  unaffected  unless  the 

additive  was  injected  immediately  in  front  of  the  antenna.  Good  results  were  also  obtained 
with  Teflon  products. 

3.6  EFFECTS  OF  GRADIENTS 

As  noted  above,  the  boundary  layer  surrounding  a  slender  cone  vehicle  in  hypersonic  flight 
has  strong  gradients  of  density,  temperature,  velocity,  and  electron  density.  In  addition, 
the  entire  flow  field  is  highly  nonuniform,  with  large  changes  in  properties  in  passing  from 
one  region  to  another.  Thus  the  importance  of  these  nonuniformities  must  be  evaluated 

before  the  breakdown  analysis  can  be  simplified. 

? 
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Fante  has  analyzed  the  situation  of  nonuniform  media  in  terms  of  two  problems.  In  the 
first  problem  he  treated  the  effect  of  an  ambipolar  diffusion  layer  and  an  outer  free  diffusion 
layer  covering  an  aperture  antenna.  The  results  of  this  calculation  show  that  even  a  small 
ambipolar  diffusion  layer  can  drastically  reduce  the  breakdown  field  strength.  The  second 
problem  c<  lsidered  is  that  of  a  hemispherical  region  in  which  the  field,  and  hence  the 
ionization  frequency,  is  zero  outside  of  an  inner  hemisphere.  Although  both  of  these  problems 
illustrate  the  importance  of  considering  nonuniform  conditions,  they  were  both  done  using 
locally  homogeneous  media  in  the  analysis.  Nevertheless,  Fante  points  out  how  to  solve  the 
more  general  problem  in  terms  of  a  generalized  electron  continuity  equation, 

If  "  **  (7Dn-nV)  +  (vrva)n  (3.6-1) 

where  n  is  the  electron  density,  D  its  the  diffusion  coefficient,  V  is  the  velocity  vector,  v.  is 
the  ioni  ation  frc  "lency,  and  V&  is  the  attachment  frequency. 


m 


In  general,  it  can  be  seen  from  the  effects  of  reentry  that  several  factors  contribute  to 
spatial  nonuniformity  of  different  terms  in  the  continuity  equation.  The  flow  velocity  is 
obviously  nonuniform  in  space  as  determined  by  the  flow  field  properties.  The  diffusion 
coefficient  will  have  gradients  due  to  particle  number  density  gradients  and  ionization 
gradients  (transit m  to  ambipolar  diffusion  in  high  electron  density  regions).  If  the  effective 
field  is  included  in  the  formula  for  the  diffusion  coefficient,  gradients  will  also  be  due  to 
field  nonuniformity.  The  frequencies  of  ionization  and  attachment  have  spatial  gradients 
contributed  by  variations  in  particle  number  density,  field  intensity,  and  also  air  tempera¬ 
ture.  Profiles  of  chemical  additives  which  are  distributed  in  the  flow  field  would  conceivably 
contribute  to  gradients  in  both  the  diffusion  coefficient  (molecular  weight  effect)  and  the 
ionization  and  attachment  frequencies.  Finally,  the  ionization  gradients  would  affect  pulse 
breakdown  in  terms  of  the  spatial  variation  of  initial  electron  density,  although  this  effect 
is  probably  not  very  important. 
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SECTION  4 

METHODS  OP  ANALYSIS  OP  BREAKDOWN 

Generally  speaking  there  are  two  methods  of  analysis  whloh  oan  be  used  to  predlot 
antenna  breakdown:  the  ktnetie  theory  method  and  the  phenomenal  theory,  in  the 
e«c  theoty  approach  an  attempt  Is  made  to  solve  the  Boltzmam.  equation  for  the 
electron  distribution  function,  using  data  on  collision  cross  section. 

*“”“**“'  meth0d  -  breakdown  has  been  widely  used  by  various 

p27n  °  ‘“‘I  breakd°Wn  ‘he0Iy-  ^  —  *  a  collection  of  extra- 

Lc  71T  °"  b88t  “  °f  reaS°MbIe  phenomeI10Iogical  ideas  to  existing 

basic^.  This  method  starts  with  foe  electron  continuity  equation,  which  can  be 

erivet,  from  foe  Boltzmann  equation.  The  various  terms  ln  this  equation  are  then 

Z  “  7Ple  °' th'  ^  WrtoblM  -  *  and  experimental 

data  are  used  to  determine  the  arbitrary  constants  in  these  equations. 

;rr whloh  approaoh  18  best  ^ for  predic“ng  the  ^  °<  ***** 

antenna  breakdown  is  very  important.  Some  subsidy  issues  need  to  be  setJbefore 
this  question  oan  be  answered.  The  following  ore  Important  questions: 

tte  JZ^rX^CfoXo'SSfo  , T  breakd0M  WI“  make 

b'  w?rk“s°ut  the  ciaasicai  lhe0,y — 

Mos*  Investlgators  of  reentry  effects  have  **  to  analyse  the  situation  ty  extending  the 

^—logical  *“*  N°  — «■-*  ba.  yet  been  encounter^  exceptor 

rrr  a  t*  d“cuity  ,s  “  •• "  °f  -  «*■*  »e  mo 

great  to  use  this  approach:  however,  the  same  uifficuliy  would  have  a  much  worse  effect 
on  the  practicality  of  foe  kinetic  theory  method 
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The  reason  the  kinetic  theory  approach  is  not  favored  classically  is  that  closed  form 
solutions  have  been  obtained  only  in  very  simple  cases  where  the  electron  collision 
frequency  is  a  simple  function  of  velocity.  However,  Carleton  and  Megill66  have 
successfully  coded  the  Boltzmann  equation  for  solution  by  a  digital  computer  for  the 
case  of  a  weakly  ionized  air  without  ionization,  attachment,  or  recombination  reactio 
Thus  the  kinetic  theory  is  a  useful  approach  when  collision  cross-section  data  are 


•  dnu  Ajenanaer  or  Aerospace  Corporation  have  adapted  the  program  of 

Carleton  and  Megill  to  the  solution  of  the  Boltzmann  equation  in  the  breakdown  problem. 
The  program  includes  the  effects  on  the  electron  distribution  function  of  elastic  collisions 
with  nitrogen  and  oxygen  molecules  and  inelastic  collisions  including  rotational,  vibra¬ 
tional,  and  electronic  excitation.  It  requires  30  seconds  to  1  minute  on  a  CPC  computer, 
per  case.  The  solution  of  the  diffusion  equation  by  a  variational  method  is  an  important 
ps  rt  of  the  calculation.  If  the  effects  of  the  effective  field  concept  are  included  in  the 
uncertainly,  the  predicted  breakdown  conditions  are  within  the  scatter  of  MacDonald's 
data  at  all  pressures,  although  the  predictions  do  not  agree  with  generally  accepted 
phenomenological  theory  at  low  pressures  (below  10  torr).  All  work  so  far  has  been 
applicable  to  cold  air  breakdown.  It  remains  to  be  seen  whether  this  program  will  give 
good  predictions  of  re-entry  effects. 


SECTION  5 


CONCLUSION 


In  view  of  the  apparent  role  of  molecular  excitation  in  the  effect  of  high  gas  temperature  on 
antenna  breakdown,  the  problem  of  predicting  the  excitation  levels  for  various  flow  field 
constituents  becomes  very  important.  Not  only  would  the  nonequilibrium  vibrational  excitation 
of  flow  field  gases  be  an  important  part  of  the  prediction  of  breakdown,  but  also  the  inter¬ 
pretation  of  experimental  data  on  high  temperature  breakdown  would  hinge  on  a  specification 
of  the  vibrational  state  of  the  test  gas.  Thus  the  solution  of  the  problem  of  vibrational 
excitation  is  important  in  deriving  theoretical  models  from  experimental  data  as  well  as  in 
applying  such  models  to  reentry  breakdown  predictions,  allhough  dissociation  may  be  niore 
important  than  vibration  in  terms  of  temperature  effect  on  breakdown. 


Based  on  an  extensive,  but 


as  vet  incomplete,  literature  survey  on  (he  effects  of  vibrational 
nonequilibrium  in  shock  layer  applications,  the  following  provisional  conclusions  may  be 
drawn: 


An  analysis  on  a  microscopic  scale  which  would  provide  the  nonequilibrium 
distributions  of  vibrational  energy  level  populations  during  the  simultaneous 
relaxation  of  vibration  and  dissociation  is,  at  present,  not  possible  for  the 
air  system. 


r  _ ... 


On  the  condition  that  the  vibrational  energy  levels  are  populated  according  ,o 
a  continuous  sequence  of  Boltzmann  distributions,  a  macroscopic  analysis  which 
offers  an  approximate  solution  to  the  coupled  vibration-dissociation-recombination 
relaxation  problem  is  possible  for  the  air  system  and  is,  in  fact,  aval  *ble. 
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